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Abstract. 

Motions of stars in close binary systems with a conservative mass exchange are ex- 
amined. It is shown that Paczynski-Huang model widely used now for obtaining the 
semi-major axis variation of a relative stars orbit is incorrect, because it brings about 
large mistakes. A new model suitable for elliptical orbits of stars is proposed. Both of 
reactive and attractive forces between stars and a substance of the flowing jet are taken 
into account. A possibility of a mass exchange at presence of accretion disk is considered. 
PACS numbers: 97.80.Fk, 97.10.Gz 



INTRODUCTION 



The research of stars motions in close binary systems began in the 60-s' of the last 
century in a cycle of works by Kruszevski [1], Hadjidemetriou [2], Piotrovski [3], Huang 
[4], Paczynski [5] and others. As a result for a case of conservative mass exchange the 
simplified dependence was received for semi-major axis a of a circular relative stars orbit 
from constant mass increase of an accepting star M 2 as 

This dependence, which we shall call for by Paczynski-Huang name and up to this 
day it is used in all researches dedicated to close binary systems with a conservative mass 
exchange when Mi + M 2 = M = const. 

The derivation of formula (1) is based on the assumption that equations of motion for 
stars with variable masses admit the angular momentum integral 

J = M1R4 x Vi + M 2 R 2 x V 2 = const, (2) 

where Mj,Rj, Vj, (i = 1,2) — the mass, radius-vector and velocity of stars motion. By 
the most an assumption is made that the close binary stars form a closed mechanical sys- 
tem which admits integrals of momentum and angular momentum. Differential equations 
of motion thus are usually not written out. 
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However the assumption on existence of integral of momentum is erroneous. We shall 
be convinced of it, proceeding from most general view of the Mestschersky equations for 
a two-body problem with variable masses [6] : 

rfVi „M l M 2 n d\ 2 „M l M 2 

M '^r = G ^ R + ^ M '-5T = - G -sr lR + *. < 3 > 

where R = R 2 — Ri, and Qi and Q 2 are jet forces acting on stars S\ and S 2 . 

Let's consider three opportunities of reception of integral of momentum (j2D from equa- 
tions (EJ. 

1. At first we shall take into account of closed mechanical system, described by the 
problem of Mestschersky-Levy-Civita [7], differential equations of which turn out from 
([3]), if jet forces put equal to: 

Q 1 = -M X V X , Q 2 = -M 2 V 2 , (4) 

where M x = -M 2 = const, M 2 > 0. 

This model is usually used in astronomy at study of motions of mutually attracted 
pair stars, taking place in dust cloud with account of jet forces, arising owing to sticking 
substances dust of cloud on stars. 

The Mestschersky-Levy-Civita problem represents a closed mechanical system, there- 
fore there are integrals of momentum and of angular momentum (121) . However to use 
model Mestschersky-Levy-Civita for study on motions of stars in close binary system 
with a conservative exchange of mass is not admitted, so as the jet forces (jlj) in close 
binary system do not exist. The true jet forces have other directions and others absolute 
values: on a donor-star the jet force acts, directed in the opposite side to radius R and 
is equal to product of velocity of sound on Mi, and on accreted-star acts the jet force, 
directed on tangent to a trajectory of relative motion of flowing particles at the moment 
of their hit on a surface of a star. These jet forces sharply differ from values (JH) in the 
Mestschersky-Levy-Civita problem. 

For this reason the model of Mestschersky-Levy-Civita is not suitable for study of 
stars motion in close binary system. 

2. Any other problems of two bodies with variable masses, admitting integral of 
angular momentum, does not exist. But it is possible to consider an opportunity of 
the approximate reception of integral (j2J) from equations ([3]). For this purpose we shall 
consider, that the jet forces acting on stars, are small values in comparison with forces of 
mutual attraction between stars, and consequently they can be neglected 

Qi =0, Q 2 = 0. (5) 

Then we shall receive a well known Gylden-Mestschersky model differential equations 
of which are 

dV, M X M 2 dV 2 Mi M 2 

Ml lu = G ^^ R ' M2 lu = - G ^^ R ' (6) 
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in which integral of the angular momentum ([2]) in absolute motion does not exist, but a 
quasi-integral does exist 

J' = M1R1 x Vi + M 2 R 2 x V 2 - J (M1R1 x Vi + M 2 R 2 x V 2 ) dt = const. (7) 

If in this quasi-integral one neglect small members, containing as multiplier the velocity 
of mass change Mi and M 2 , we shall receive integral ([2]) , but thus a problem of Gylden- 
Mestschersky will be transformed to a two-body problem with constant masses, with the 
help of which it is impossible to study a transfer of substance between stars. 

Nevertheless the conservative (M = const) problem of Gylden-Mestschersky represents 
the certain interest for study conservative exchange in close binary system, so as it is 
possible with the help of it to characterize approximately motions of stars at presence of 
accretion disk. The equations of relative motion in the Gylden-Mestschersky problem are 

~dt ~ "1^ R (8) 

and at M — const exactly coincide with the equations of relative motions of the classical 
two-body problem with constant masses [8] , therefore for relative motion of bodies in the 
problem of Gylden-Mestschersky there is a strict angular momentum integral 



GMa(l - e 2 ) = const, (9) 

and all Kepler's elements of an orbit in this problem are constants, including a = const. 

Thus, the exclusion from consideration all of jet forces, i.e. the transition to a conser- 
vative problem of Gylden-Mestschersky, results to invariance of the semi-major axis of a 
relative orbit of stars, but does not result in the formula of Paczynski- Huang. 

3. We shall consider now the third opportunity of a "formal conclusion" the rules ([T]), 
which is usually used in the literature. Proceeding from the lack of external forces acting 
on a system, and invariance of its complete mass, one asserts, that the system is closed 
and, hence, admits the existence of integral (TS]). 

Such statement is valid for stars with constant masses, but is not valid for stars with 
variable masses. The point is that for systems of bodies with variable masses to be closed 
except of absence of external forces and the invariance of complete mass is needed it 
is necessary a presence some additional "internal" jet forces, which make system to be 
closed. 

Really, according to the second law of Newtons mechanics the velocity speeds of change 
of momentum of stars are determined by equalities: 

rf(MiVi) r M 1 M 1 d(M 2 V 2 ) Mi M 2 

where F\ and F 2 are external forces, acting on stars. 
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Let's rewrite these equations as 



Mi 



i 



G 



M X M 2 
R 3 



R + Q1 + F1, M 2 



dV 2 



-G 



M 1 M 2 
R 3 



R+Q 2 + F 2 , 



(11) 



dt 



dt 



where Qi = —M{V\ and Q 2 = — M 2 V 2 are "internal" jet forces, active on stars Si and 
S 2 . 



From there it is visible, that integral of momentum M{V i + M 2 V 2 = const and 
together with it those of angular moment (j2J), will exist at absence of external forces, 
Fi = F 2 = 0, but it is obviously necessary a preservation of "internal" jet forces Qi = 
-MiVi and Q 2 = -M 2 V 2 . 

Therefore for existence of integral (|2J) it is necessary not only the absence of external 
forces and invariance of complete mass of the system, but also the presence of "internal" 
jet forces of a kind fll]) too is necessary. These requirements are carried out only for the 
problem of Mestschersky-Levy-Civita, considered in item 1. For a problem on motion of 
stars in close binary system with a conservative exchange of mass these requirements are 
not carried out, so as the jet forces have another vector values which are distinct from 
"internal" jet forces, namely, Qi = — MiWi and Q 2 = — Af 2 W 2 , where Wi ^ Vi and 
W 2 V 2 are relative velocities of the outflow (inflow) of mass on stars. Therefore the 
integral of angular momentum does not here exist. 

Thus, the assumption on existence of integral of angular momentum (J2J) is infaithfull, 
i.e. those about a problem on stars motion in close binary system with a conservative 
exchange of mass being closed, is erroneous. In a consequence of that the model of 
Paczynski- Huang is incorrect, and it cannot be used for study of motions in close binary 
stars. As show results of numerical integration of equations of motion with taking into 
account of true jet forces, acting on a star, the use of Paczynski- Huang model results 
in significant mistakes in definition of the semi-major axis, down to an opposite sign of 
derivative a. 

For this reason any conventional correct model, determining motions of close binary 
stars, now does not exist. By the aim of purpose our work is the such model creation. 
For circular motions of stars such a model was offered in the work [9], in which except 
for forces of mutual attraction of stars both the true jet forces, and forces of attraction 
on stars from jet stream of substance are taken into account also. The definition of stars 
motions is carried out with the help of numerical integration. In the present work a similar 
research, but under assumption, that the orbits of stars are elliptic will be carried out. 



For definition of motion of flowing particles of masses we shall use the restricted ellip- 
tical three-body problem, the equations of which in rotating and pulsating, barycentric 



THE RESTRICTED ELLIPTICAL THREE-BODY PROBLEM 



4 



system of coordinates of Schepner (Scheibner, Petr, Nechvil, Rein) x,y look like [10] 

d 2 x ^ dy dU 

dv 2 dv ^ dx ' / 19 \ 
d?y 2 dx_dU_ ^> 

dv 2 dv ^ dy ' 

where the x-axis is always directed on a accreted-star S 2 , v is true anomaly of stars, 
p = 1/(1 + ecosf) is the dimensionless distance of a star S 2 comparatively S±. Through 
U the Jacobi function in the restricted elliptic three-body problem is designated 

x 2 + y 2 3 f l-m m\ p 2 2 

U = — - — + p + — )+— (3 + m -m), 13 

2 V r\ r 2 J 2 

and r\ and r 2 are distances from a flowing particle of a jet accordingly up to the centers 
of masses of the first and second stars 



T\ = \J (x + pm) 2 + y 2 , r 2 = \J (x + pm — p) 2 + y 2 , (14) 

where 1 — m — M\jM and m = M 2 /M are relative masses of stars, M = Mi + M 2 , 
p = a(l— e 2 ) is the focal parameter of an orbit, a and e are semi- major axis and eccentricity 
accordingly 

The stream of substance from a donor-star occurs through a vicinity of internal Euler 
libration point L±, xl with which is conclude in limits — pm < xl < p — pm and is 
determined numerically as a root of the nonlinear equation dU/dx = at y = 0, i.e. of 
equation 

, / 1 — m m \ , . 

x-p 3 [^^=+ - =0. 15 

V v x + P m v x + P m ~ ' v 7 

Roche lobes in the planar elliptic restricted three-body problem are determined with 
the help of equations of curves of minimal energy [11] 

1 -p 2 (3 + m 2 -m) = C(l + ecosw), (16) 

r± r 2 J 

where C is Jacobi's constant. 

Let's consider, that the accepting star has the form of sphere 

p 2 [(x + pm - p) 2 + y 2 } = P 2 , (17) 

the radius P of which in stream process of substance changes according to the dependence 



P = Po ?/— • (18) 

The stream of substance begins with a donor-star after achievement by flowing parti- 
cles a level of energy more than that of in the lipration point L\. As it has been shown in 



5 



[11], the stream of substance through the vicinity of a point L\ has a pulsating character 
and occurs in the vicinity of apoastr of orbits. The velocity Vq of the stream of particles 
masses from a star S\ is always directed along the x-axis to the star S 2 and has the value 



IGM 



V Q (v) = VaoViiv), V x {v) = m^——Vl + 2ecosv + e 2 J (19) 

where Voo — 0.03 is the coefficient, established from observations. The apoastr vicinity, in 
which occurs the stream of mass, is determined on a level of energy of particles (Jacobi's 
constant C, see [11]) by range of true anomaly tt — v a < v < n + v a , where < v a < n. 
Therefore the mean velocity V c of stream of mass from the star Si during one orbital 
period of stars is possible to determine with the help of formula 



K+ r Va fGM 1 *Y a 

/ V dv = V 00 mJ / v / l + 2ecosw + e 2 dv. (20) 

2v a J V P 2v a J 



All stream lines of a jet, having scatter both on coordinates and on time during one 
"emission" of substance for an orbital period, we shall approximate by one trajectory, 
outgoing from singular point L\. The numerical integration of equations ffl2|) is carried 
out with such initial conditions: 

v = ir, v a = -, x(vq) = x L , x'(v ) = -p y(vo) = 0, y'{v ) = 0, 
m = mo, fa = const, P = Pq, 

where "prime" signifies differentiation on true anomaly, (. . .)' = d(. . .)/dv. 

The numerical integration is carried out on an interval vq < v < vq + r, where the 
value of true anomaly v + r is determined by the moment of hit of a particle on the 
second star surface in the point x 2 = x(v + r),y 2 = y(v + r). If such a moment does 
not exist, the process of stream of substance occurs with forming of an accretion disk. 

As a result of numerical integration there are components of outflow velocity of Wi 
from the star Si and the inflow mass velocity W 2 on the star S 2 : 

W lx = V c , W ly = 0, W 2x = x'(v + T )v, W 2y = y'(v + t)v. (22) 

Jet forces acting on stars Si and S 2 are considered to be applicable to their centers of 
masses and are determined with these formulas 

Qi = mM{-V c , 0}, Q 2 = mM{x'(v + T)v,y , (v + r)v}. (23) 

With the help of numerical integration also the mass of established stream S3 and 
coordinates X3, y% of its center of masses are determined as: 

VO+T VO+T 

M 3 = —mM, X 3 = - f x(v)dv, Y 3 = - I y(v)dv, (24) 

V c T J T J 

vo vo 
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where 



2tt 2tt 



2 



If., If GM .„ n2 , GM ( e 

— I v dv = — I \ — — ( 1 + e cos v ) dv = \ — — 1 H 

2nJ 2tt J V P 3 V P V 2 





is the mean angular velocity of orbital motion. 

With the help of formulas (T2"5j) and (1211) the jet forces, the mass of jet and coordinates 
of its center of masses are determined. These values enter in the right-hand sides of the 
differential equations of relative motion of stars. 

DIFFERENTIAL EQUATIONS OF MOTION OF STARS 

As initial differential equations of motion of stars in an inertial coordinate system we 

shall consider system of equations: 

dVx MjMo MiM 3 
Mi— = G — R + Ch + G— — r 13 , 

dV 2 M M, M 2 M 3 1 J 

M 2 — = -G^^R + Q 2 + 



where r 13 = y (x 3 + pm) 2 + y 3 , r 23 = y (#3 + P m — p) 2 + 2/1 distances between centers 
of masses of stars and that of the jet. The equation of a star 5 2 relative S\ motion from 
there turns out to be as 

£ = ™ R+ *_* +GM ,fqi_!j»), (26 ) 

dt i? 3 Af 2 Mi 3 V r 23 r is/ 

where V = V 2 — Vi is the velocity of a star 5* 2 relative to S\ motion. 

Differential equations in osculating elements for semi-major axis of a relative orbit a 
and its eccentricity e are represented as [10]: 

— = , [e sin v S + (1 + ecosf )X1, 

dt v /GMa(l - e 2 ) (27) 



<ie / a(l — e 2 ) 



smvb + I cosu -I ) 1 

1 + e cos v 



dt V GM 

Perturbing accelerations S and T are happened to be known after integration of equa- 
tions O: 



c . /W^x , V c \ m fx 3 +pm-p x 3 +pm 

S — m\ h + GMt— ; ; 

m \ — mJ v c V r 23 rj 3 



(2? 



23 ' 31 



Choosing as the independent variable the true anomaly v of stars, the equations ([2 
can be rewritten as follows 

da 2a 2 



^=[esinvS + (1 + ecos v)T], 

dv y^ GMa(l- e 2 ) 

de 1 a(l — e 2 ) f . / cos t> + e \ " 

-r- = tI/ Sinto+ COSH + - i 

dv GM [ y 1 + ecosw 
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where 



' 1 + e cost>) + " 



p 3 v ' V GM 



1 + e cos f 



After averaging on the true anomaly v, equations of stars motion result in the form 



2n 2 2 

/ — . \e sin v S + (1 + e cos v)T]dv, 



da 1 r 2a 

L V / — — ( 30 ) 

de 1 / 1 a{l -e 2 ) \ . / cosv + e\ 1 

— = — / -\/ — — smv b + cosv H 1 dv. 

dv 2n J v\ GM [ V 1 + ecosW J 

o 

Besides if instead of v we choose the new independent variable m with the help of 

d Tfi d 

the formula — = : — , then equations of motion for averaged elements a and e will 

dv i) c dm 

happened to be as 

da v r r 2a 



- — = — — [ — . fe sin vS + (1 + e cos v)T]dv, 

dm 2vrm7 v ^GMa{l - e 2 ) 



de v c [ 1 \a{\ — e 2 ) 



(31) 



dm 2nm J v V CM 



o 



sm do + [ cos v H ) 1 

1 + e cos f 



dv. 



After calculation of definite integrals in right-hand sites of equations (131]) we shall 
receive the final form of differential equations for determination of secular perturbations 
of semi-major axis and eccentricity of the stars relative orbit: 



(im (1 — e 2 ) 2 m 
rfe 2 + g 2 e W 2w 3 



(32) 



rfm 25(1 - e 2 ) 1 + Vl - e 2 m 2 



ry 3 r 3 



where the designation W 2y = y'(v + r) is accepted. 

If to put Wi y = 0, i.e. not to take into account jet forces, the equations (I3"2l are 
simplified and accept the form 

da 



e 2 ry 3 r 3 , 



— = 2d 3 VT 

dm ^ 

— = -ea 2 Vl - e 2 ry 3 r 3 , 
dm 2 

which admin the existence of such first integral 

a 3 e 4 = const = a 3 , (34) 

with initial conditions do = a(mo) and eo = e(mo). 

Further for simplicity the overlined index "tilde" above values a and e will be omitted. 
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NUMERICAL RESULTS 



The relative orbit of stars was determined by means of numerical integration of dif- 
ferential equations ( |32|) and ( 133]) . The integration was carried out on an interval by 
independent variable m from 0.1 up to 0.8. 

In fig. 1 are given dependencies of the semi-major axis a and eccentricity e of relative 
orbit of a accreted-star from value q = m/(l — m) for various values of initial radius stars 
Po at initial value of eccentricity eo = 0.1. From this figure it is visible, that the change of 
an orbit strongly depends on Pq. At small values Pq the flowing jet of substance gets on 
the surface of a star S2 with the large velocity and also creates enough the large jet force 
accelerating motion of a star, while for large Pq, on the contrary, the jet force is directed 
in the opposite side and consequently drags the motion of the star. 

The influence of initial value of eccentricity eo on elements of star S2 is represented 
in fig. 2 for fixed value P/a = 0.075. As it is seen stars motions strongly depend on the 
initial value of eccentricity eo- Under small eo and Pq secular changes of the semi- major 
axes can reach significant values. 

In fig. 3 on the plane (q,P/a) the curves are represented describing a close binary 
system. The curve 1 determines border P max /a — 0.49 \/m, after which the binary forms 
contact system, when the sizes of accreted-star reach the singular Euler point. The solid 
curves 2-5 for various eo determine border for P m i n /a below of which at close binary 
system an accreted disk is formed. Between curves Pmax/a and P m i n /a the half-divided 
phase of close binary system settles down. On dot and dash curves 6-9 for various eo the 
first derivative of the semi-major axis turns into zero. Above these curves a close binary 
system extends, and downwords it is compressed. The model of Paczynski-Huang on this 
figure is represented by a vertical straight line q — 1. According to the formula ( [1]) at 
q < 1 a binary is compressed, and at q > 1 it extends, that sharply differs from results 
which has been received by means of numerical integration. 

In fig. 4 on the same plane (q,P/a) curves of change are represented for value P/a 
depending on q for various values Pq. The curve 1 crosses the border of contact binary 
system formation. The further account for this curve is not carried out. The curves 6 
and 7 cross a border of accreted disk formation. If one considers, that the accreted disk 
is formed instantly, it is possible to continue these curves further, using the differential 
equations (1531) . i.e. believing, that any jet stream does not create jet force on a star S 2 , 
giving back energy on rotary motion for particles of accreted disk. Continuations of curves 
6-7 on the figure is shown by dashed lines. 

In fig. 5 the comparison of change of the semi-major axis and eccentricity of an orbit 
of stars is carried out at presence of accreted disk and at its absence for various initial 
values of eccentricity. As it is visible from figure changes a and e at presence of accreted 
disk are not rather great, that is explained by absence of a jet force and small forces of 
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attraction on stars by the stream. Therefore for, approximate estimations of evolution of 
close stars orbit with presence of accreted disk, taking into account small sizes of a jet 
and, hence, small right hand sides of equations fl33|) . it is possible to consider, that 

a = const, e = const, (35) 

i.e. instead of Paczynski- Huang model one can use that of Gylden-Mestschersky. 

Dependencies between the semi-major axis a end eccentricity e along trajectories of 
stars motion are represented in fig. 6. For the accreted disk these dependencies turn out 
from the first integral flMj) . 

It is difficult to carry out comparison of the present work results with those of other 
authors, so as jet forces, created by a stream, actually nobody takes into account (see, for 
example, [11]), but in the present work the jet forces play the basic role. 

It is doubtless, that all received results require a further revision and their conformity 
to observations. 

CONCLUSION 

During the last half-century for definition a relative orbit of close binary stars the 
incorrect Paczynski- Huang model was used, that was marked in the paper [8]. However 
up to this day in works, dedicated to close binary systems this model continues to be used. 
Therefore in the present work once again, but by means of other methods, the inaccuracy 
of this model has been shown. 

For determination of relative motion of stars in close binary system in the present work 
the numerical integration of equations of motion is used with taking into account of jet 
forces and forces of attraction of stars by the flowing jet. The above calculations of elliptic 
orbits of close binary stars show that the eccentricity of orbits can change strongly enough. 
The influence of jet force on orbital evolution of stars can be various. If the accepting star 
has smaller mass and little sizes, and the accreted disk is absent, then the eccentricity is 
increased strongly, reaching values 0.5-0.8. Probably just the such influence of jet force 
explains rather large number of close double stars with great orbital eccentricity. If the 
accepting star has the large mass and the great size, the jet force creates the braking 
effect, and the eccentricity decreases down to zero. 

It is shown that for approximate determination of orbital evolution of close binary 
system with a generated accreted disk instead of the Paczynski-Huang model it is neces- 
sary to use the model of Gylden-Mestschersky, according to which all Keplerien elements 
remain constant, in particular, a = const and e = const. 

It is doubtless, that obtained results can be specified, taking into account other per- 
turbation factors and making new assumptions based on observations. 

The present work is maintained by the Russian fund of fundamental researches (grant 
08-02-00398). 
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(a) 



(b) 



Figure 1: The diagrams of functions: (a) — a(q) and (b) — e(q) at a = 1 and e = 0.1 
for various P . Notations: 1 - - P = 0.06, 2 — P = 0.073, 3 — P = 0.075, 4 - 
P = 0.0758, 5 — P = 0.0765, 6 — P = 0.1. 
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Figure 3: Characteristic curves on a plane yq, — J. Notations: 1 — the border of forma- 
tion of the contact system, solid curves 2, 3, 4, 5 — the border of accreted disk formations 
for values eo = 0, eo = 0.1, eo = 0.3, eo = 0.5 accordingly, doted and dashed curves 6, 7, 8, 
9 — the border of change of a sign of derivative a' for the same values eo accordingly, verti- 
cal direct line 10 — the border of change sign by derivative a' under the Paczynski-Huang 
formula. 
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Figure 4: Dependence of the ratio — from q for different trajectories of motion of stars 
at a = 1 and e = 0.1. Notations: \ — P = 0.1, 2 — P = 0.0878, 3 — P = 0.08, 4 
- P = 0.0778, 5 — P = 0.0765, 6 — P = 0.0758, 7 — P = 0.075, 8 — the border 
of contact system, 9 — the border of accreted disk. Dashed curves — continuation the 
appropriate curves after formation of accreted disk, doted end dashed curve — the border 
of change of sign of derivative a'. 
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1 2 ' 3 ' 4 1 2 3 4 



(a) (b) 

Figure 5: The diagrams of functions: (a) — a(q) and (b) — e(q) for a — 1 and P = 
0.0758. Notations: — e = 0, 1 — e = 0.1, 2 — e = 0.3, 3 — e = 0.5. Solid curves 
correspond to absence of accreted disk, and dot end dash curve — to its presence. 




(a) Without of accreted disk at do = 1 and cq = (b) With accreted disk. Notations: 1 — eo = 0.1, 

0.5. Notations: 1 — P = 0.0758, 2 — P = 2 — e = 0.3, 3 — e = 0.5. 
0.075, 3 — P = 0.073. 

Figure 6: Dependencies between a and e along trajectories. 
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